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ABSTRACT 

The  initial-value  problem  is  studied  for  evolution  equations  in  Hilbert 
space  of  the  general  form 


A(u)  +  8(u)  a  f 

where  A  and  B  are  maximal  monotone  operators.  Existence  of  a  solution  is 
proved  when  A  is  a  subgradient  and  either  A  is  strongly-monotone  or  8  is 
coercive?  existence  is  established  also  in  the  case  where  A  is  strongly- 
monotone  and  8  is  subgradient.  Uniqueness  is  proved  when  one  of  A  or  8 
is  continuous  self-adjoint  and  the  sum  is  strictly-monotone?  examples  of  non¬ 
uniqueness  are  given.  Applications  are  indicated  for  various  classes  of 
degenerate  nonlinear  partial  differential  equations  or  systems  of  mixed 
elliptic-parabolic-pseudoparabolic  types  and  problems  with  non-local 
nonlinearity. 
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SIGNIFICANCE  AND  EXPLANATION 


A  variety  of  models  of  physical  phenomena  can  be  written  in  a  unified  way 
as 

—  Au  +  8u  3  f 

dt 

where  A  and  8  are  monotone  operators  between  Hilbert  spaces.  Therefore  it 
seems  convenient  to  study  questions  of  existence  and  uniqueness  for  abstract 
equations  as  the  above  in  view  of  their  range  of  applications. 

He  will  divide  the  applications  in  two  categories:  problems  with  local 
non-linearities  and  problems  with  global  non-linearities.  The  former  occur 
when  at  some  point  of  the  region  where  the  process  occur,  the  function 
describing  the  process  either  jumps  or  degenerates.  Typical  are  the  Stefan 
problem,  diffusion  in  porous  media,  diffusion  in  partially  saturated  porous 
media.  They  include  also  certain  kinds  of  diffusion  in  a  medium  with  a 
singularity  due  to  a  crack  or  fissure. 

Problems  with  global  non-linearities  are  typical  of  processes  where  a 
threshold  (which  separates  the  different  aspects  of  the  phenomenon,  say  for 
example  elliptic-parabolic)  is  given  through  the  global  measure  of  the  energy 
(or  the  variation  of  the  energy). 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


IMPLICIT  DEGENERATE  EVOLUTION  EQUATIONS  AND  APPLICATIONS 


Emmanuels  DiBenedetto 


and  R.  E.  Showalter 


1.  Introduction. 

Let  A  and  B  ba  maximal  monotone  operators  from  a  Hilbert  space  V  to  its  dual 

V*.  Such  operators  are  in  general  multi-valued  andtheir  basic  properties  will  be  recalled 

below.  He  shall  consider  initial-value  problems  of  the  form 

(1.1)  7-  IA( a))  +8( a)  3  t  ,  u(0 )  j  v 

at  0 

2  *  • 

where  f  e  L  (0,TjV  )  and  Vg  e  V  are  given.  It  is  assumed  throughout  our  work 

that  A  is  a  compact  operator  from  v  to  V*.  In  applications  to  partial  differential 
equations  this  assumption  limits  the  order  of  the  operator  A  to  be  strictly  lower  than 
that  of  B.  Both  operators  will  be  required  to  satisfy  boundedness  conditions  and  one  or 
the  other  is  assumed  to  be  a  subgradient. 

The  objective  of  this  work  is  to  prove  existence  of  a  solution  of  (1.1)  when  A 
and  8  are  possibly  degenerate.  Observe  that  we  must  in  general  assume  some  condition  of 
coercivity  on  the  pair  of  operators.  To  see  this,  we  note  that  if  one  of  them  is 
identically  zero  then  (1.1)  is  equivalent  to  a  one-parameter  family  of  "stationary" 
problems  of  the  form  M(u(t))  3  F(t)  where  M  is  maximal  monotone.  But  if  M  is,  e.g., 
a  subgradient  in  a  space  of  finite  dimension,  it  is  surjective  only  if  it  is  coercive. 

Thus  it  is  appropriate  to  assume  at  least  one  of  A  or  8  is  coercive.  In  accord  with 
this  remark  our  work  will  proceed  as  follows.  First  we  replace  A  by  the  coercive  operator 
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A  +  eR,  where  e  >  0  and  R  :  V  •»  V  is  the  Riesz  isomorphism  determined  by  the  scalar 
product  on  V,  and  we  solve  the  initial-value  problem  for  the  "regularized"  equation 
(1.2)  7-  (A  +  eR)(u  >  +  B(u  )  »  f  . 

Here  we  may  take  e  -  1  with  no  loss  of  generality  and  we  make  no  coercivity  assumptions 
on  either  A  or  5.  Next  we  assume  B  is  coercive  and  let  e  +  0+  in  order  to  recover 
(1.1)  with  (possibly)  degenerate  A.  Since  R  is  of  the  same  order  as  B  this 
regularization  is  analogous  to  the  Yoshida  approximation.  The  operator  A  is  assumed  to 
be  a  subgradient  in  the  above.  Finally,  we  show  the  initial  value  problem  can  be  solved 
for  (1.2)  when  8  (but  not  necessrily  A)  is  a  subgradient. 

He  mention  some  related  work  on  equations  of  the  form  in  (1.1).  The  theory  of  such 
implicit  evolution  equations  divides  historically  into  three  cases.  The  first  and 
certainly  the  easiest  is  where  8  o  A  1  is  Llpschitz  or  monotone  in  some  space  [6,  23]. 
The  second  is  that  one  of  the  operators  is  (linear)  self-adjoint,  and  this  case  includes 
the  majority  of  the  applications  to  problems  where  singular  or  degenerate  behavior  arises 
due  to  spatial  coefficients  or  geometry  (2,  25).  These  situations  are  described  in  the 
book  [9]  to  which  we  refer  for  details  and  a  very  extensive  bibliography.  The  third  case 
is  that  wherein  both  operators  are  possibly  nonlinear.  This  considerably  more  difficult 
case  has  been  investigated  by  Grange  and  Mignot  [12]  and  more  recently  by  Barbu  [4].  In 
both  of  these  studies  a  compactness  assumption  similar  to  ours  is  made.  Our  boundedness 
assumptions  are  more  restrictive  than  those  in  the  papers  above,  but  they  assume  f  is 
smooth  and  that  both  operators  are  subgradients.  By  not  requiring  that  S  be  a 
subgradient  in  (1.1)  we  obtain  a  significantly  larger  class  of  applications  to  partial 
differential  equations,  especially  to  systems. 

Our  work  is  organized  as  follows.  In  Section  2  we  recall  certain  information  on 
maximal  monotone  operators  and  then  state  our  results  on  the  existence  of  solutions  of  the 
initial-value  problems  (1.1)  and  for  (1.2).  The  proofs  are  given  in  Sections  3  and  4. 
Section  5  contains  elementary  examples  of  how  non-uniqueness  occurs,  and  we  show  there  that 
uniqueness  holds  in  the  situation  where  one  of  the  operators  is  self-adjoint.  Section  6  is 
concerned  with  the  structure  and  construction  of  maximal  monotone  operators  between  Hilbert 


.1 


-2- 


t 


2.  Preliminaries  and  Main  Results. 

We  begin  by  reviewing  information  on  maximal  monotone  operators.  Refer  to  [1,  3,  11] 
for  additional  related  material  and  proofs.  Then  we  shall  state  our  existence  theorems  for 
the  Cauchy  problem  (1.1). 

Let  V  be  real  Hilbert  space  and  A  a  subset  of  the  product  V  x  v.  We  regard  A 
as  a  function  from  V  to  2V,  the  set  of  subsets  of  V,  or  as  a  multi-valued  mapping  or 
operator  from  V  into  V}  thus,  f  *  A(u)  means  [u.f]  «  A.  We  define  the  domain 
D( A)  ■  (u  (  Vi  to  non-empty},  range  R  (A)  -  U  {aui  u  *  V)  and  inverse 

g 

A-1(u)  ■  {v  «  V:  u  e  A{ v) }  of  A  as  indicated.  The  operator  A  is  monotone  if 

(f1  -  fj,  u1  -  “2*v  *  0  whenever  *  *  for  J  -  1,2.  This  is  equivalent  to 

(I  +  XA)-1  being  a  contraction  for  every  X  >  0.  We  call  A  maximal  monotone  if  it  is 

maximal  in  the  sense  of  inclusion  of  graphs.  Then  we  have  a  monotone  A  maximal  monotone 

if  and  only  if  R  (I  +  XA)  -  V  for  some  (hence,  all)  X  >  0.  If  A  is  maximal  monotone 
9 

we  can  define  its  resolvent  3  (I  +  XA)  ,  a  contraction  defined  on  all  V,  and  its 

Yoahlda  approximation  A^  -  X-1(I  -  J^),  a  monotone  Lipschitx  function  defined  on  all  V. 

For  ut  V  we  have  A^(u)  e  A(J^(u)).  We  denote  weak  convergence  of  xn  to  x  by 

x  -  x. 
n 

Lemma  2.1.  Let  A  be  maximal  monotone,  [x_.v_]  e  A  for  n  >  1,  x  -  x.  v  *  v  and 
' "  n  n  n  n 

11m  inf(yn,xn)v  «  (y,x)y.  Then  (x.yl  «  A.  If  in  addition  lim  sup(yn,xn)v  «  (y,x)y, 

then  (y  ,x  )„  ♦  (y,x)...  We  observe  that  A  Induces  on  L2(0,T»V)  a  maximal  monotone 
n  n  v  v 

operator  (denoted  also  by  A)  defined  by  v  t  A(u)  if  and  only  if 
v(t)  e  A(u(t))  for  a.e.  tr  (0,T). 

A  special  class  of  maximal  monotone  operators  arises  as  follows.  If  :V  *  (-*,“]  is 
a  proper,  convex  and  lower  semlcontlnuous  function  we  define  the  subgradlent  !vC  V  »  V  by 
3*(x)  •  {*  e  V:  *(y)  -  *>(x)  >  (*,y-x)  for  all  y  r  v}  . 


The  operator  is  maximal  monotone.  Furthermore  it  is  useful  to  consider  the  convex 


The  following  are  equivalent:  z  e  3^(x),  x  e  3  0*(z),  and  i^(x)  +  **(z)  «  (x,z)v;  thus 
^3  is  the  inverse  of  SV’.  we  mention  the  following  chain  rule,  tl].  Let  H'tO.TjV) 
denote  the  space  of  absolutely  continuous  V-valued  functions  on  fO,T]  whose  derivatives 
belong  to  L2(0,T»V). 

Lemma  2.2.  If  u  «  H^O.TiV),  v  e  L2(0,TjV)  and  [u(t),v(t)]  e  3^  for  a.e.  t  e  [0,T], 
then  the  function  t  +  ^  (u(  t) )  is  absolutely  continuous  on  [0 ,T]  and 

<p(u(t>)  -  (w,u'(t^)  ,  all  w  *  3i^(u(t))  , 

for  a.e.  t  e  [0,T] . 

There  is  a  version  of  monotone  operator  from  V  to  its  dual  space  V*  which  is 

•  # 

equivalent  to  the  above  through  the  Riesz  map  R  :  V  ♦  V  .  Thus,  AC  v  x  V  is  monotone 
if  and  only  if  A  =  R  1  o  A  is  monotone  in  V  x  V  and  maximal  monotone  if  and  only  if 
Rq( R  ♦  A)  •  V*  in  addition.  We  shall  use  these  two  equivalent  notions  interchangeably. 

Our  applications  to  partial  differential  equations  will  lead  to  operators  on  V  x  V*.  Also 
the  subgradient  is  naturally  constructed  in  the  W  -  W*  duality  of  a  Banach  (or 
topological  vector)  apace  W  .  Finally  we  cite  the  following  chain  rule. 

Lemma  2.3.  Let  V  and  W  be  locally  convex  spaces  with  duals  V*  and  w*.  Let 

*  *  * 

A  :  V  *  W  be  continuous  and  linear  with  dual  A  :W  ♦  V  .  If  •f  :  w  ♦  is 

proper,  convex  and  lower  semicontinuous  then  so  also  is  o  A  :  V  ♦  ,  and  if  is 

continuous  at  some  point  of  R  (A)  we  have  [11] 

9 

* 

3(*  o  A)  «  A  o  o  A  , 

Our  results  on  the  existence  of  solutions  of  the  Cauchy  problem  (1.1)  are  stated  as 
follows. 

Theorem  1.  Let  W  be  a  reflexive  Banach  space  and  V  a  Hilbert  space  which  is  dense  and 

imbedded  compactly  in  W.  Denote  the  injection  by  i  :  V  ♦  W  and  the  dual  (restriction) 

*  *  * 

operator  by  i  :  w  ♦  V  .  Assume  the  following: 

[A^]  The  real-valued  V  is  proper,  convex  and  lower  semicontinuous  on  w,  continuous  at 
some  point  of  V,  and  o  i  :  v  +  W  is  bounded. 


[B^]  The  operator  6  :  V  ♦  v  is  maximal  monotone  and  bounded. 

Define  A  s  i  o  X  o  i,  Then  for  each  given  f  £  L2(0,T>V*)  and  tu0,v0)£  A  there 
exists  a  triple  u  £  H’fO.T.-V),  v  £  H’tO.TjV*),  and  w  e  L2(0,T;V*)  such  that 
<2.1. a)  Jj-  (Ru(t)  +  v(t) )  +  w(t)  =  f(t)  , 

( 2 . 1 ,b)  v(t)  £  A(u(t)>,  w(t)  f  8(u(t)),  a.e.  t£  [0,T]  , 

(2.1.c)  Ru(0)  +  v(0)  =  Ruq  +  vQ  . 

Theorem  2.  In  addition  to  the  above,  assume 

[A2]  3*  o  i  :  L2 (0 ,T; V )  -  L2(0,T;W*)  is  bounded 

,2  2  * 

[B2 )  t>  :  L  (0,T;V)  *  L  (0,TfV  )  is  bounded  and  coercive,  i.e.. 


II  ull 


lim 

L‘(0,T,V) 
tu,  v]£  B 


/  Tv(t)(u(t))dt 
0 

innr - 

L  <0,T;V) 


Then  for  each  given  f  e  L2(0,T;V  )  and  Vg  e  Rg(A)  there  exists  a  triple 
u  £  L2(0,T;V),  v  £  H1 ( 0 ,T jV*) ,  w  £  L2(0,TiV*)  such  that 
(2. 2. a)  v(t)  +  w(t)  “  f(t>  , 

( 2 . 2.b)  v(t)  £  A(u(t)),  w(t)  £  B(u(t)),  a.e.  t  E  [0,T]  , 

(2.2.c)  v(0)  -  Vp. 

Remarks.  From  Lemma  2.3  it  follows  that  A  »  3(^lv>  where  *|v  -  if  o  i  is  the 
restriction  of  ■f  to  V.  Since  A  :  V  ♦  V  is  bounded  it  follows  D(A)  «  V,  hence, 

VC  D(3s5)  C  dom(^)  C  w 

and  y  is  continuous  on  the  space  V.  Also,  since  , (0)  <  "  we  may  assume  with  no  loss 

of  generality  that  r(0)  <  0  and  thus  y  (z)  >  0  for  all  z  £  V. 

*  *  *  * 

From  the  compactness  of  i  :  W  *  V  it  follows  A:  V  +  V  is  compact,  i.e.,  maps 
bounded  sets  into  relatively  compact  sets. 

Since  B  is  bounded  and  maximal  monotone  we  have  D(8)  «  V.  It  is  important  for  our 


applications  that  we  have  made  no  assumptions  which  directly  relate  A  and  8 
Specifically,  we  do  not  compare  A(x)  and  8(x)  in  angle  or  in  norm. 


Finally,  we  give  a  variation  on  Theorem  1  in  which  only  the  second  operator  8  is  a 


subgradient.  The  compactness  assumption  on  A  is  retained. 

Theorem  3.  Let  the  spaces  v  and  w  be  given  as  before.  Assume  the  following: 

*  *  * 

(Aj]  The  operator  A  :  V  ♦  V  is  maximal  monotone  with  R^(  A)  C  w  and  A:  V  +  W 

bounded. 

tB j]  The  real-valued  i(i  is  proper,  convex  and  lower  semicontinuous  on  V  and 
8  =  3i|»  :  V  ♦  V  is  bounded. 

Then  for  given  f  e  L2  (0,T;V*)  and  (u0,v0Jf  A  there  exists  a  triple  u  e  H1  (0,T;V), 
v  e  H^O.T.-V*),  and  w  f  L2(0,T;V*)  satisfying  (2.1). 


3.  Proofs  of  Theorem  1  and  Theorem  3 


These  proofs  are  very  similar;  let  us  consider  first  the  Theorem  1.  He  formulate 
(2.1)  in  the  space  V.  Set  a  «  R  2  o  A#  b  *  8  '  o  8,  etc,  and  consider  the  equivalent 
equation 

(3.1. a)  (u(t)  +  v( t) )  +  w(t>  -  f(t>  , 

(3.1  .b)  v(t)  £  A( u( t) ) ,  w(t)  e  B( u( t) ) ,  a.e.  t  e  (0,T)  . 

Let  X  >  0  and  consider  the  approximation  of  (3.1)  by 

(3. 2. a)  (u^(t)  +  (v^(t) )  +  Bx(ux(t))  -  f(t)  , 

(3.2 .b)  vx(t)  e  A(ux(t)),  t  £  [0 ,T)  . 

Since  (I  +  A)-1  and  bx  are  both  Lipschitz  continuous  from  V  to  V,  (3.2)  has  a 
unique  absolutely  continuous  solution  ux  with  ux(0)  +  vx(0)  "  u0  +  v0‘  sl-nce 
(I  +  A)-1  is  a  function,  we  have  ux(0)  “  u0  and  vx(0)  "  vo* 

we  derive  a  priori  estimates  on  ux .  Take  the  scalar  product  in  V  of  (3. 2. a)  with 
ux<t)  and  note 

(V*(t),ux(t))v  -  /(vx) 

by  Lemma  2.2  where  is  the  conjugate  of  ^lv  in  V.  Integrating  the  resulting 

identity  gives 

Vfcllux(t)lly  +  <t  (vx<t)  ^V'jllugll2  +  *  <v0> 
t 

♦  /  (If  f  ( s)II  v  +  llBx(0)llv)llux(8)flvds,  0  <  t  <  T  . 

Since  {B^ (0 ) }  is  bounded  by  the  fact  that  0  £  D(B) ,  *  >0  and  f  r  L2(0,T;V),  we  have 

proved  the  first  part  of  the  following. 

Lemma  3.1.  The  following  are  bounded  independent  of  X  >  0: 


II  UX" 


(a) 


L  (0,TjV) 


I  Rvxii  „  .  , 

L  (0 ,Tj W  ) 


1  Jx^ux>" 


L  (0,T;V) 


"VV11 


L  (0 ,T; V) 


(b) 


X  L2(0,T;V) 


11  v*n  2 

A  L*(0,T;V) 
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Proof:  The  second  and  third  terms  of  (a)  are  bounded  because  the  operators  A  :  V  ♦  w 
and  J  5  (I  ♦  Xb)-1  :  V  ♦  V  are  bounded.  Since  £  B*Jx<ux"  and  B  is  bounded' 

the  last  term  in  (a)  is  bounded. 

To  obtain  (b)  we  take  the  scalar  product  of  (3. 2. a)  by  u^(t),  note  that 
(v^(t),  u^(t) >  >0  by  (3.2.b)  and  monotonicity  of  A,  and  thereby  obtain 

lu'(t)«^  <  (»ftt)lv  +  IBx(ux(t)llv)«u^(t)«v  , 
so  we  bound  the  first  term  in  (b).  The  second  follows  from  (3. 2. a). 

Note  that  we  have  {Rvx}  bounded  in  L2(0,T;W*)  and  {Rvj^}  bounded  in 
L2(0,T;V*).  Since  W*  is  compact  in  V*  it  follows  from  [17,  p.  58]  that  {Rvx)  is 
(strongly)  relatively  compact  in  L  (0,T;V  ).  From  this  observation  and  Lemma  3.1  it 
follows  we  may  pass  to  a  subsequence,  again  denoted  by  ux»vx,  for  which  we  have 


(3. 3. a) 


( 3 . 3 .b) 


(3.3.C) 


UX  *  U'  W  *  UX  *  U>  ' 

2 

vx  ♦  v  (strongly),  vj^  v'  in  L  (0,T;V)  , 
ux<t)  *  u(t)  and  vx(t)  *  v(t),  all  t  t  [0,T] 


Since  ux  -  Jx<ux>  =  Xbx(ux>  ♦  0  there  follows 


( 3 . 3 • d) 


Jx(ux)  -  u  in  L  ( 0 ,T : V)  . 


It  remains  to  show  that  u,  v,  w  satisfy  (3.1)  and  the  initial  condition.  First  we 
use  (3. 3. a)  and  (3.3.b)  and  Lemma  2.1  to  obtain  v  e  A(u).  Next  take  the  scalar  product  of 
(3. 2. a)  with  any  x  e  V  and  integrate  to  get 


(ux<t)  +  vx(t),x)v  +  /  (Bx(ux( s) ) ,x)vds  =  /  (f(s),x)vds  +  (uQ  +  vQ,x)^ 


Takina  the  limit  as  X  -»  0  gives  (since  x  is  arbitrary) 


u(t)  +  v(t)  +  /  (w  -  fids  =  ^  +  Vg ,  0  <  t  <  T 


From  this  identity  we  obtain  (3.1. a)  and  u(P)  +  v(0)  =  Up  +  Vp ;  since  v(0)  f  A(u(0)) 
and  (I  +  A)  1  is  a  function  we  have  u(P)  =  up.  In  order  to  show  w  e  B(u),  and  thereby 
finish  the  proof  of  Theorem  1,  it  suffices  hv  Lemma  .’,1  to  show 


lim  sup  (B.(u.),J  (u.  )  -  u)  , 

VO  *  X  L  lO,T>V) 


<  0  . 


We  note  further  that 


(W'  Jx<uxn  "  (W'  W  -  V  +  'VV'V 


-X(Bx(Ux),  W*  +  <VV'V 


so  it  suffices  to  show 


(3.4) 


lim  sup  (B.(u.),u.  -  u) 

X*0  *  A  L2(0,T,V) 


<  0  . 


By  (3. 2. a)  it  follows  (3.4)  is  equivalent  to 


(3.5) 


lim  inf  (u*  +  v>  u.  -  u) 

X-*0  L2(0,T»V> 


>  0 


Define  t|i(x)  -  ^  #x*2  +  ^(x),  x  f  V  so  that  3iJ>  -  1  +  3*.  From  (3.2.b)  and  Lenrnia 
2.2  we  obtain 

( u^  (t)  +  V^(t),ux(t))v  -  I|«*(ux(t)  ♦  vx(t))  , 

and  by  integrating  there  follows 

(uj  +  v’,  u  )  «  iMu.(T)  +  v.  (T) )  -  **(u„  ♦  v„>  . 

A  A  A  L2(0,T;V)  XX  00 


Similarly  we  have  from  (3.1. a) 
(u*  +  v’ ,u) 


»  i(i  (u(T)  +  v(T) )  -  ^  (u  ♦  v  ) 
L  (0,T;V)  0  0 


By  (3.3.c)  and  weak  lower  semicontinuity  of  tji  we  have 

*  * 

<|i  ( u(T )  +  v(T ) )  <  lim  inf  *  <u.(T)  +  v.(T))  , 
X*0  A  A 


and  our  preceding  calculations  show  this  is  equivalent  to  (3.5). 

Remark  3,1.  From  Lemma  2.1  we  find  that 

<B.(u  ),  JX(UX,)  ■>  *  (w'u)  ■> 

L2(0,T»V)  L  (0,TiV) 


If  we  also  have  R  (or  6)  strongly  monotone  then  we  can  take  the  limit  in  the  estimate 


»  clJ. (u.  >  -  u»2 

L  (0,T:V) 


'****?• 


<VV  -  w,  Jx(ux)  -  u) 

L  (0,T;V) 

to  conclude  {Jx(uxJ  and  {ux}  converge  strongly  to  u  in  L^(0,T;V). 

Remark  3.2.  It  is  clear  that  we  actually  have  v(t)  e  A(u(t))  for  every  t  f  [0,T). 

The  proof  of  Theorem  3  closely  follows  the  preceding  pattern.  That  is,  formulate 

(2.1)  as  the  equivalent  initial  value  problem  for  (3.1)  and  approximate  this  by  (3.2)  with 

u, (0)  +  v  (0)  »  u.  +  v.  for  each  X  >  0. 

X  X  0  0 

To  derive  a  priori  hounds  we  take  the  scalar  product  of  (3. 2. a)  with  ux(t)  and 
integrate  to  obtain 


(3.6) 


T  T 

/  lu'.v  +  t|lx(Ux(T))  <  +  /  (f(t).U'(t))vdt  . 


Here  ^  is  the  Yoshida  approximat ion  of  We  may  assume  iJj  is  non-negative  and  the 

same  holds  for  ,  so  we  have  the  first  part  of  the  following. 

Lemma  5.1.  The  following  are  bounded  independent  of  X  >  0: 

(a)  *V-  ,  2  -  Wv A»  .  .  . 

L  (0 ,T;V )  L  ( 0  ,  T ;  V )  L  (0,T;W  ) 


(b) 


«J.(u,)ll  ,  llB,(u,)ll 

L*(0,T;V)  l"(0,T;V) 


V  -  -  llvi"  2 

L  (0  ,T;  V  )  L  (0  ,T;V ) 


Proof:  The  bound  on  the  first  two  terms  in  (a)  follow  from  (3.6)  and  the  remaining  terms 
in  (a)  are  bounded  by  [Aj]  and  [Bjl.  Next  take  the  scalar  product  of  (3. 2. a)  with  vx(t) 
and  obtain  (b)  as  was  done  in  Lemma  3.1. 

We  may  pass  to  a  subsequence  satisfying  (3.3)  and  we  obtain  as  in  Theorem  1  the  triple 
u,  v,  w  satisfying  the  equation  (3.1. a)  and  initial  condition  and  v(t)  e  Au(t), 
t  e  (0 ,T] •  It  remains  to  show  w  f  B(u)  and  this  is  equivalent  to  showing  (c.f.  (3.6)) 
lim  inf  (u!  +  vj,u.)  > 

X-*n  X  A  A  Lz(0(T;V) 
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(3.7) 


>  (v',u) 


L2(0,T,-V) 


Since  uj^  e  L  (0,TrV)  we  may  integrate  by  parts  to  compute 
(3. 8. a) 


(UX  +  vA'V  2  "  -  'Vv> 

L  <0,T;V)  A  V  0  V 


-<v.,u!)  ?  +  (V.  (T),u.  (T))„  -  (v.,u_ 

x  x  L2(0.T,V)  X  X  V  0  0 


and  similarly,  since  u'  e  L2(0,T»V), 


(o'  ♦  v',  u* ) 


(3.8.b) 


L  <0,T»V) 
-  ( V,u*  ) 


t^)(«u(T)|2  -  luQl2) 


2  +  <v(T),u(T))  -  (V  ,u  ) 

I<  (0,T|V)  V  0  0  V 


Finally  we  observe  that  (3.7)  follows  immediately  from  (3.3)  and  (3.8). 


4.  Proof  of  Theorem  2 


Choose  un  e  a'^Vq).  For  each  X  >  0  let  u^,v^  e  hVo.TjV),  w^  t  L2(0,T»V) 
satisfy 

(4.1. a)  Xu^(t>  +  vx(t)  ♦  wx(t)  “  f<t)  ' 

(4.1.b)  v^(t)  *  A<u^(t)),  w^(t)  f  B(u^(t)),  a.e.  t  e  [0,T]  , 

(4.1.c)  Xu^(O)  +  v^(0)  «  XuQ  +  v0  . 

The  problem  (4.1)  has  such  a  solution  by  Theorem  1,  and  our  plan  is  to  show  that  we  may 
take  the  limit  as  X  *  0  in  (4.1)  to  obtain  a  solution  u,w  e  L2(0,TiV),  v  e  h'(0,T?V)  of 
(4. 2. a)  v* (t)  +  w(t)  -  f(t)  , 

(4.2. b)  v(t)  e  A(u(t)>,  w(t)  t  B(u(t))»  a.e.  t  e  (0,T)  , 

(4.2.C)  v(0)  ”  vq  • 

With  our  notation  A  *  R  1  o  A,  etc.,  (4.2)  is  equivalent  to  (2.2). 

We  proceed  to  derive  a  priori  estimates.  Consider  first  the  initial  condition.  Since 
(XI  +  A)  '  is  a  function  it  follows  from  (4.1.c)  that 


(4.3)  UX(0)  *  u0'  VX<0)  "  v0'  X  >  0  . 

Lemma  4.1.  The  following  are  bounded  independent  of  X  >  0  •• 

(a)  llu.<  ,  X^llu.H  „  , 

A  L*( 0,T;V)  L  ( 0 ,T;V ) 


We  drop  the  second  (non-negative)  term  in  (4.4)  and  note  by  monotonicity  of  B  that 
(w^,  u^)v  >  (£»u^)v  for  some  5  f  B(0).  Thus  (4.4)  gives 


13- 


+  c 


I  <wrVv<af, 


L  (0 ,T(V)  "  L  (0 ,T)V ) 


and  the  coercivity  of  B  implies  the  boundedness  of  the  first  term  in  (a).  The  second  now 
follows  from  (4.4)  and  now  part  (b)  follows  from  our  assumptions  A1  and  . 

Lemma  4.2.  The  following  are  bounded  independent  of  X  >  0: 


"VX"  2 

L  <0,T:V) 


IXuM 

L  (0,T;V) 


Proof:  Take  the  scalar  product  of  (4.1. a)  with  vjjt).  Since  (u^(t),  vj^(t))  >0  by 

monotonicity  of  A  we  obtain 

#vj|(t)»y  <  (»f(t)#v  +  «wx(t)llv)IVx(t)lv 

from  which  the  first  bound  is  immediate.  To  obtain  the  second  we  take  the  scalar-product 
of  (4.1. a)  with  u^(t)  and  <3f°P  the  non-negative  term  (ux (t) ,vx<t) This  gives 

X #ux (t) <  (»f(t)«v  +  5wx(t)»v)Iu> (t)«v 
and  hence  the  desired  bound. 

We  have  now  shown  that  {Rvx}  is  bounded  in  L2(0,T;W*)  and  that  {Rv’}  is  bounded 
in  L2(0,TiV*).  Since  W*  is  compact  In  V*  it  follows  that  (Rv,^)  X®  strongly  compact 

•y  * 

in  L  (0,T?V  ).  From  this  observation.  Lemma  4.1,  and  Lemma  4.2  it  follows  we  may  pass  to 
a  subsequence  (which  we  denote  again  by  {ux},  {vx},  {wx})  for  which  in  L2(0,TtV)  we 


v,  *  v,  v.'  -  v' 


Note  that  Xux  ♦  0  and  it  follows  Xux  -  0  by  standard  arguments.  Furthermore  we  may 
assume  vx(t)  *  v(t)  in  v  for  a13-  e  tO,T]  by  equi-continuity  of  {vx},  and 

similarly  Xux(t)  *  0  in  V  for  all  t  e  [0,TJ. 

It  remains  to  show  that  the  triple  u,  v,  w  obtained  above  constitutes  a  solution  of 

> 

(4.2).  Let  x  f  V,  take  the  scalar  product  of  x  with  (4.1. a)  and  Integrate  to  obtain 


t  t 

(Xu^(t)  +  v^(t),x)v  +  /  (w^(s),x)vds  -  /  ( f  ( s) ,  x)  v<is  +  (XuQ  +  vQ  ,x>v 


Since  weak  convergence  in  LZ(0,T;V)  implies  weak  convergence  in  L^(0,tiV).  it  follows  by 
letting  X  ■*  0  that 


t  t 

(v(t),x)  +  /  (w(s),x)  da  ■  /  (f(s),x)  ds  +  v  ,  x  e  V,  t  e  [0,T)  . 

"ft  "ft  VO 


That  is, 


t  t 

v(t)  +  /  w(s)ds  »  /  f(s)ds  +  v  ,  a  .e.  t  e  tO,T]  , 

0  0 

and  this  implies  (4. 2. a)  and  (4.2.c).  From  Lemma  2.1  there  follows  v  e  A(u)  so  it 
remains  only  to  establish  we  B(u).  For  this  it  suffices  by  Lemma  2.1  to  show 

(4.4)  lim  sup  (w.,u  )  <  (w,u)  . 

X+0  LZ (0 ,T;V)  LZ(0,TjV) 

In  order  to  prove  (4.4)  we  first  note  by  (4.1. a)  and  (4. 2. a)  that  it  is  equivalent  to 

(4.5)  lim  inf  (Xu*  +  v',u  )  >  (v',u) 

X-»0  A  A  LZ (0 ,T; V  >  L  (0,T;V> 

-1  *  * 

Since  ux(t)  f  A  ^X^*'  =  ^X^^  on  [0,T]  ,  where  •f  is  the  conjugate  of 

•f  |^,  we  obtain  from  Lemma  2.2 

(Xu^(t)  +  vjk(t),ux(t))v  -  ((X/2)lux(t)1y  +  <#>*(  vx(  t) ) } 

so  we  integrate  and  obtain 

( XuJ  +  v!,u.)  ,  =  <X/2)lu.(T)!2  +  /(V.  (T>>  -  (X/2XU.I2-  /(vn) 

A  A  A  LZ(0,T,V)  X  V  X  n  V  0 

>  v'*(vx(T))  -  (X/2)lu0H2  -  /(vQ)  . 
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J 


Similarly  we  compute 


(V  ,u) 


-  if  ( v(T) )  -  f  (v.)  . 
L  <0,T(V)  U 


Since  {v^}  are  equi-uniformlycontinuous  we  have  v^(t)  ♦  v(t)  at  every  t  e  tO,T],  so 

the  lower  semicontinuity  of  <f  *  gives 

lim  inf  *  ( v.  (T ) )  >  ^  ( v(T) )  . 

\*0  A 


In  view  of  the  preceding  computations  this  is  exactly  (4.5), 


5.  Remarks  on  Uniqueness. 

We  first  present  an  example  which  shows  that  gross  non-uniqueness  of  solutions  of 
(1.1)  can  occur,  even  if  both  operators  are  strongly  monotone  subgradients.  Moreover  the 
non-unioueness  occurs  in  each  term  of  the  triple,  u,  v,  w,  not  just  in  the  latter  two 
terms  selected,  respectively,  from  A(u)  and  B(u).  Next  we  shall  show  that  uniqueness 
does  hold  for  (1.1)  when  at  least  one  of  the  operators  is  continuous,  linear  and  symmetric 
and  the  sum  of  the  operators  is  strictly  monotone.  Our  last  example  shows  that  symmetry  of 
the  linear  operator  is  essential. 

Example  1.  Let  V  «  w  ■  R,  the  space  of  real  numbers,  and  define 

A(s)  -  B(s)  -  s  ♦  H(S-1)  , 

(  1  ,  r  >  1 

where  H(r)  -  <  (0,1)  ,  r  »  0  , 

[  0  ,  r  <  1 

denotes  the  Heaviside  function  and  f  =  0.  Consider  the  initial-value  problem  (1.1)  which 
takes  the  form 


r  v* (t)  +  wit)  -  0,  v(0)  -  2 

(5.1) 

i  v(t>  -  u(t)  f  H(u( t )  -  1),  w(t)  -  u(t)  e  H(u(t) 

-  1) 

» 

Let  g 

he  any  maximal  monotone  graph  or  continuous  function  from  R  to  R 

such 

that 

g(s)  = 

s  for  s  <  (1,2)  and  g(s>  C  [1,2)  for  s  e  (1,2).  Then,  if  v 

is  a 

solution 

of 

(5.2) 

v'(t)  +  gfvft))  “  0  ,  t  >  0,  v(0)  -  2  , 

it  follows  that  with  u(t>  =  A  '(v(t))  and  w(t)  5  -v'(t)  we  have  a  solution  of  (5.1). 

This  procedure  yields  an  abundance  of  solutions. 

We  display  some  special  cases  of  the  above.  Pick  c  €  [V^,1]  and  define  gc  to  bo 

that  maximal  monotone  graph  such  that  g  <t)  ■  {c  1 } ,  t  e  (1,2),  and  a  (t)  “  { t }  , 

c  o 

t  /  (1,2).  The  corresponding  solution  vc  of  (5.2)  is  given  by 

v  (t)  =  2  -  t/c,  0  <  t  <  c,  v  <t)  -  eC“r , 
c  c 
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t  >  c  . 


I 


J 


With  the  two  functions  uc  and  wc  given  by 

u  (t)  «  1,  w  (t)  ”  1/c  for  0  <  t  <  c 
c  c 

u  (t)  »  w  (t)  »  eC_t,  t  »  c  , 
c  c 

this  provides  a  continuum  of  solutions  of  (5.1). 

We  can  give  the  following  elementary  sufficient  conditions  for  uniqueness  to  hold  for 

(1.1)  or,  equivalently,  for  (4.2). 

Theorem  4.  Let  A  and  S  be  monotone  operators  on  a  Hilbert  space  V.  Suppose  A  ♦  B 
is  strictly  monotone  and  that  one  of  A  or  B  is  continuous,  linear  and  symmetric.  Then 
for  each  function  f  :  [0,T]  ♦  V  and  vQ  e  v  there  is  at  most  one  solution  u,  v,  w  of 

(4.2) . 

Proof:  Suppose  A  is  continuous,  linear  and  symmetric.  For  j  »  1,2  let  u  ^ ,  v^  Wj  be 
a  solution  of  (4.2).  Take  the  scalar  product  of  the  difference  of  (4. 2. a)  with  u^  -  Uj 
to  obtain 

(Ifc  )  (A(Uj  (t)  -  u2<t)),u1(t)  -  u2(t))v  ♦  (wi(t)  -  Wjltl.u^t)  -  u2(t))v  »  0 
Integrating  this  identity  and  using  (4.2.c)  gives 

t 

(Vi  )(A(u1(t)  -  u  (t)),u  <t)  -  u  (t))y  +  /  (w1  -  w2,ui  -  u2>vds  -  0,  0  <  t  < 

0 


and  this  implies 

Au  ^  ( t )  =  Au  2  ( t ) ,  (w^t)  -  w2(t),  u^t)  -  u2(t))v  =  0  a.e.  t  e  (0,T). 
Since  A  +  B  is  strictly  monotone  we  have  u^ft)  =  Uj(t),  hence 
v^lt)  -  Au^t)  =  Au2(t)  =  v2(t),  and,  by  (4.1. a),  w^(t)  »  w2(t)  a.e.  on  (0,T). 

Suppose  now  B  is  continuous  linear  and  symmetric.  Starting  with  two  solutions  as 
above  we  integrate  the  corresponding  equations  (4. 2. a)  to  obtain 


t 

v . ( t)  +  B(9  . ( t ) )  -  v.  +  /  f ,  j 

1  3  On 


,2 
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IV 


(5.3) 


vwaswwf* . ,v«^tV'-  ^»fev-tqfe^^p»irg^i  n  n  i  n  >  1 1  s  i  |  nmmiiii  Iiinfiyill 


t 

where  0  (t)  =  /  u..  Takinq  the  difference  of  (5.3)  for  j  *  1,2,  then  scalar-product 

3  0  3 

with  0J  -  6^  and  integrating  gives  us 
t 

<5-4>  /  (v,  -  v2,  e*  -  o^)v  +  0/2  )(B(ot(ti  -  o2(tn,  e^t)  -  e2<t))v  =  o  . 

Since  v^(t)  f  A(9’(t))  a.e.,  each  term  is  non-negative.  It  follows  that 

B(8|(t)  -  Sjft))  ”  0  on  !0,T]  and  thus  from  (5.4)  that 

(v^t)  -  v2(t),  u^t)  -  u2(t)>v  =  0  a.e.  t  f  (0,T)  , 

so  the  desired  results  follows  by  strict  monotonicity  of  A  +  R. 

Finally  we  cite  an  example  to  show  that  the  symmetry  condition  cannot  be  eliminated 

\ 

from  Theorem  4. 

Example  2.  Let  H^(0,1)  be  the  Sobolev  space  of  those  absolutely  continuous  functions  on 
the  interval  (0,1)  whose  first  derivatives  belong  to  L2(0,1);  set 
V  «  {v  e  H  (0,1)  :  v(1)  =0)  and  note  that  V  C  L2 (0,1)  C  V*.  Define  A  :  V  *  V*  by 

A(v)  =  -v* .  Clearly  A  is  linear  and  we  have 

1  , 

( v )  (  v )  =  -/  v'v  ds  *  C/2  )  I  v  ( 0  )  |  >  n 

0 


so  A  is  monotone.  Let  $  he  aiven  by 


and  define 


f 

\ 


r/2 , 


r2/2. 


r  <  0  or  r  >  1  , 
n  <  r  <  1  , 


1 

6(u)(v)  =  j  ?(u'(s))v*(s)ds, 
0 


u ,  v  <. 


V  . 


It  is  easy  to  check  that  t>  is  a  strictly  monotone  suhqra-Uent  on  V. 


(S.5) 


Consider  the  Cauchy  problem 

(A(u) )  +  8(u)  -  0,  Au(0)  -  -1 

with  the  above  operators.  A  solution  u  of  (5. 5)  is  a  weak  solution  of 
boundary-value  problem 

(5. 6. a)  <-u  >  -  (8(u  )>  -  0,  0  <  x  <  1,  0  <  t  , 

X  t  XX 

<5.6.b)  ux(0,t)  -  u( 1  ft)  -  0  , 

(5.6.C)  -ux<x,0)  »  -1 

where  the  subscripts  denote  partial  derivatives.  Consider  the  following 


(1), 

u  <x,t) 


U<2,<x.t> 


f  <X2  +  t2)/2t  -1,  0  <  x  <  t  <  1  , 
[  x  -  1  ,  0  <  t  <  x  <  1  , 
t/2  -  1,  0  <  x  <  t/2  <  V2  . 


I  x  -  1,  0  <  t/2  <  x  <  1 ,  t  <  1  . 

It  is  a  straightforward  computation  to  check  that  both  u^  ^  and  2  ^ 
(5.6),  hence,  both  are  solutions  of  (5.5).  Note  that  the  only  condition 
met  in  this  example  is  the  symmetry  of  A.  It  shows  also  that  8  being 
not  a  satisfactory  substitute  for  8  to  be  continuous  and  self-adjoint. 


the  initial- 


two  functions: 


are  solutions  of 
of  Theorem  4  not 
a  subgradient  is 


6.  Construction  of  Differential  operators. 

We  have  been  discussing  evolution  equations  which  contain  a  pair  of  nonlinear 
operators  from  a  Hilbert  space  V  to  its  dual  V*.  In  our  applications  the  generalised 
solutions  obtained  in  our  theorems  may  satisfy  natural  or  variational  boundary  conditions 
(e.g.,  of  Neumann  type)  which  are  implicit  in  the  functional  identity 
(6.1)  dt  A(u(t))  +  f<fc> 

in  V*.  Such  boundary  conditions  are  classically  recovered  by  Green's  formula  so  we  shall 
describe  an  appropriate  extension  of  this  formula  which  requires  a  minimum  of  regularity  of 
the  generalized  solution.  The  objective  is  to  resolve  each  term  in  (6.1)  into  two  parts,  a 
differential  operator  in  distributions  over  a  region  0,  the  formal  operator,  and  a 
constraint  on  the  boundary  T,  the  boundary  operator.  Then  we  briefly  recall  basic  facts 
on  Sobolev  spaces  and  construct  a  rather  general  nonlinear  operator  8  which  will  be  used 
in  the  next  section  to  illustrate  our  Theorems  in  some  examples  of  initial-boundary-value 
problems. 

Assume  we  are  given  a  linear  surjection  y  i  V  *  T,  called  a  "trace"  operator,  which 
is  a  strict  homomorphism  onto  its  range  T,  called  "boundary  values"  of  V.  Let  Vg  be 
the  kernel  of  y  and  note  that  the  dual  operator,  y  ( g)  »  g  o  y,  is  an  isomorphism  of 
the  dual  space  T*  onto  the  annihilator  in  V*.  Suppose  there  is  given  a  continuous 

seminorm  |*|  on  V  for  which  Vg  is  dense  in  the  seminorm  space  U  5  {V, I  - 1 } .  Then  we 

#  *  * 

naturally  identify  U  simultaneously  as  a  subspace  of  V  and  of  Vg. 

V*  * 

We  resolve  the  operator  A  :  V  *  2  into  a  formal  part  in  Vg  and  a  boundary  part 

in  T*.  For  each  u  e  D[AJ  set  A(u)»  !  F  e  A(u)},  the  set  of  restrictions  to  Vg 

0 

of  functionals  in  A(u).  Then  set  D[Agl  5  (u  f  V  :  Ag(u)  n  u  *  $}  and  define 
Aq  :  v  ♦  2U  by  Ag(u)  *  Ag(u)  n  That  is,  Ag  is  the  set  of  those  functionals  in 

Ag(u)  which  have  (unique)  continuous  extensions  in  U*  C  V*.  Now  let  u  e  D(Ag]  and 


F  e  A(u) 

with  Fg  «  F|  e  U  ,  hence, 

c 

F0 

e  Ag(u).  Then  in  Vg  we  have 

i 

o 

N 

•  * 

Y  ( g)  for  a  unique  g  f  T  , 

• 

so  we  can  define  3^(u)  c  T  to  be  the 

set  of 

all  such 

g.  Thus,  for  each  F  e  A(u) 

for 

* 

which  F  *  F|  e  U  ,  there  is  a 

0  vn 

unique 

g  e  T  for  which 


and  we  Indicate  this  by 

(6.25  A(U)  =  Aq(u)  ♦  ue  D[A01  . 

* 

Ir.  our  applications  VQ  is  a  space  of  distributions  over  £2  and  T  is  the  space  of 
boundary  values  of  the  Sobolev  space  V,  so  (6.2)  is  the  abstract  Green's  formula  for  the 
operator  A. 

In  many  examples  the  solutions  of  (6.1)  will  have  the  additional  regularity  properties 
described  below. 

\  * 

Lemma  6.1.  Let  v  e  H  (0,T;V  )  with  v(t)  e  A(u(t))  a.e.  on  (0,T],  and  set 

Vg(t)  =  v(t)lv  for  each  t  e  (0,T).  Let  Vg(t)  £  U*  and  define  g(t)  £  T*  by 

®  *  • 

v(t)  »  vQ(t)  *  Y  (9(t))  for  tf  [0,T).  If  vg(t)  f  U  a.e.  on  tO,T]  then 

g  £  H1  (0,T;T*)  and 

v'(t)  «  v^(t)  +  Y  (g’(t)),  a.e.  t  e  [0,T]  . 

The  preceding  situation  occurs,  for  example,  in  the  case  of  linear  symmetric  A  and  in 
certain  other  special  cases  (2,  9,  17,  25], 

* 

g 

Suppose  the  operator  A  is  given  as  above  and  let  a  second  operator  B  :  V  +  2  be 
given.  Resolve  it  likewise  into  two  parts, 

(6.3)  8<u)  =  8g (u)  +  y*( Jg(u) ) ,  u  £  DtBQl  . 

Let  there  be  given  fQ  e  L2(0,T?U*),  gg  e  LZ(0,T;T*),  vQ  f  [ Ag ]  and  gp  e  T*  with 

v  +  y  (g  )  e  R  [A.].  Consider  a  solution  of  the  Cauchy  problem 
0  0  g  0 

A( u( t) )  +  8 ( u( t) )  3  fg(t)  +  y*(gn(t)),  a.e.  t  £  fO,T!  , 

* 

A(u(0))  j  vQ  +  Y  (gQ)  , 
that  is,  a  triple  u,  v,  w  for  which 

f  v(t)  f  A(u(t)),  w(t)  f  8(u(t))  , 

! 

(6.4)  j  v'lt)  +  w!t)  =  fgtt)  +  Y  (gg(t)),  a.e.  t  e  (0,T)  , 

[  v(0)  =  vfl  +  Y*(g0)  • 

By  restricting  the  above  functionals  to  Vq  we  obtain 


1 - 


1 


t 

l 


I  * 
I 


(6.5) 


f 


vQ(t)  *  Aq ( u ( t ) ) >  wQ(t)  e  BQ(u(t))  , 

v^(t)  +  wQ(t)  =  f Q ( t )  in  V0,  a.e.  t  e  [0,T]  , 

V0<0)  =  vo  • 


If  Lemma  6.1  applies  then  we  obtain  wn(t)  £  u  and  t^ie  identities  (6.2)  and  (6.3)  give 


(6.6) 


r 


g.(t>  e  3  (u( t) ) ,  g  ft)  £  3_fu(t)) 

A  A  8  B 

* 

g‘(t)  +  g  (t)  =  g  (t)  in  T  ,  a.e.  t  e  [0,T]  , 

A  S  0 

gA(0)  =  g0  . 


Thus  (6.4)  implies  (6.5)  and,  in  the  situation  of  Lemma  6.1,  also  (6.6),  so  we  call  a 
solution  of  (6.4)  a  weak  solution  of  the  pair  (6.5),  (6.6).  The  first  will  give  a  partial 
differential  equation  and  the  second  yields  variational  boundary  conditions  in  our 
examples. 

Let  ft  be  a  bounded  open  set  in  Rn  which  lies  locally  on  one  side  of  its  smooth 

i  2 

boundary  T.  H  ( £2 )  is  the  space  of  functions  in  L  (SI)  for  which  each  of  the  partial 
3s?  2 

derivatives  D  ,<f  =  belongs  to  L  (ft),  1  <  j  <  n.  Letting  Dg  denote  the  identity  on 

2  i  1 

L  (S2),  we  can  express  the  norm  on  H  (ft)  by 


n  2  V 

mi  .  «  <  I  uo  Mr  )  2* 

H  (ft)  j=0  3  L  (ft) 

We  shall  let  V  be  a  closed  subspace  of  H 1  ( £1 )  containing  C^fft)  and  let  Y  :  V  ♦  L2(T) 
be  the  indicated  restriction  to  V  of  the  trace  map  (19).  We  let  T  be  the  range  of  Y 
(a  subspace  of  l/2  ( T ) )  and  denote  the  kernel  by  VQ  s  H^fft).  Since  T  is  smooth  there 
is  a  unit  outward  normal  vector  n(s)  =  [n1 ( s) , . . . ,nn( s) )  at  each  point  s  e  T.  Note  that 

as  * 

the  test  functions  Cg(ft)  are  dense  in  V0  so  the  dual  VQ  is  the  space  of  (first  order) 

distributions  on  ft.  We  refer  to  [19]  for  information  on  these  Sobolev  spaces. 

Specifically,  we  shall  use  the  trace  operator  between  Sobolev  spaces  of  fractional  order. 

* 

We  shall  construct  an  operator  8  :  V  ■>  2  which  will  occur  in  many  of  our 

examples.  For  each  integer  k,  -1  <  k  <  n,  let  there  be  given  a  continuous,  convex 

function  ib  :  R  +  R  whose  subgradient,  0,  =  ‘*'1',.'  satisfies 
k  k  K 


(6.7)  I w I  <  C(|s|  +1)  if  we  Bk<«).  s  e  R,  -1  <  k  <  n  , 

where  C  Is  some  large  constant*  Then  define  i|i  t  V  ♦  R  by 


<i(u)  “I  /  4  (D u(x))dx  +  /  (y(u(s)))ds,  u  e  V 

k-o  n  *  r 


From  the  estimates  (6.7)  it  follows  that  ij;  is  a  sum  of  continuous  convex  functions  so  we 
can  compute  its  subgradient  term-by-term.  Recall  that  the  subgradient  F  of  the  convex 
function  v  — *  /  ^  (v)dx  at  w  e  L2(!5)  is  determined  by  F(x)  e  6.(w(x)),  a.e.  x  e  15. 

n  2k 

Since  :  V  *  L  (15)  is  continuous  linear,  the  subgradient  of  the  convex  function 
v  — *  Jf  sS^fD^vJdx  at  u  e  V  is  given  by  {D*F  i  F  e  S^ID^u)  a.e.}.  See  {11,  pp.  26-28] 
and  (1,  p.  47]  for  proofs  of  these  facts.  These  observations  show  that  the  subgradient  of 

i|<  is 


(6.8) 


n  *  * 

8  (u)  -  3i>(u)  -  l  Dk6k<Dku)  *  Y  S_1(Yu),  u  e  V 


k-0 


To  be  precise,  we  have  F  e  8(u)  if  and  only  if  there  exists  e  in  L  (15), 

0  <  k  <  n,  and  f_1  e  6_1  (Yu)  in  L2(D  for  which 


F(v)  ■  /  J  f  (x)D  v(x)dx  +  /  f  ( s) v( s) ds ,  ve  v 

15  k-0  K  r 


By  restricting  the  above  to  v  f  VQ  »  HQ(15)  we  see  the  formal  part  is  the  distribution 


F|V  “  “  l  Dkfk  +  f0  e  V0  ’ 

v0  k»1  *  *  0  0 


We  denote  this  by  the  equality  (of  sets) 


(6.9) 


Bq(u)  -  -  l  DxBk<Dku>  +  8  (u) 


k»1 
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¥ 


Let's  interpret  (6.3)  with  0*  -  L2(8).  Firet,  if  D^f^  £  U*  for  1  <  k  <  n,  then  by  the 
classical  Green's  theorem  we  have  from  above 

n 

F(v)  -  F !  (v)  -  /  {  l  f  (s)n.  (s)  +  f  .(sjWtsJds,  v  e  V  . 

v0  r  k-1  k  *  _1 


Thus  u  e  0(8)  and  we  have  shown 


n 

1^  +  f_!  e  3  with  fk  e  Bk<Dku)  . 


That  is.  when  the  terms  are  as  regular  as  indicated  we  have 


(6.10) 


36(u) 


*  W,nk  4  8-1fu> 

k«1 


Furthermore,  3g  (u)  is  defined  without  these  regularity  assumptions  on  the  individual 

terms;  it  is  sufficient  to  have  F|y  e  u*.  Finally,  we  note  that  from  (6.7)  it  follows 

0 

that  6  satisfies  the  assumptions  [8,1  of  Theorem  1  and  [B21  of  Theorem  3.  It  is  also 
bounded  from  L2(0,T;V)  to  L2(0,T;V#)  and  it  will  satisfy  [B21  of  Theorem  2  if,  in 
addition  there  is  a  pair  of  numbers  K,  c  >  0  such  that 

r  2 

i|>k(s)  >  c|sl  -  K,  s  e  R,  1  <  k  <  n 
and  one  of  the  following: 

(6.11)  \  (a)  the  estimate  holds  for  k  »  0,  or 

(b)  the  estimate  holds  for  k  -  -1,  or 
_ (c)  v  f  V  and  v  ■  constant  imply  v  =  0  . 


From  (6.11)  we  can  show  that 

2 

ij<(v)  >  c,  lv«v  -  K,  ,  v «  v  , 
and  this  implies  the  coercivity  condition  in  [Bj] . 
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7.  Examples  of  Partial  Differential  Equations. 

We  shall  describe  some  examples  of  initial-boundary-value  problems  for  partial 
differential  equations  to  Illustrate  the  applications  of  our  results*  These  examples  were 
chosen  merely  to  suggest  a  variety  of  problems  that  can  be  resolved  by  our  Theorems,  and 
they  are  not  intended  to  be  best  possible  in  any  sense. 

(a).  Elliptic-Parabolic  Equations.  For  k  *  0  and  -1  ,  let  :  R  ♦  R  be  convex  and 

continuous  with  subgradient,  ,  satisfying 

!w|  <  C(  |  s  1  +1)  if  w  c  a^(s),  s  e  R  . 

Set  w  -  Hr(fl),  V2  <  r  <  1,  V  =  H 1  (fl) ,  and  note  that  V  +  w  is  compact  and 

2 

Y  :  W  ♦  L  (D  is  continuous  [19]  .  Thus  we  can  define  by 

(v)  2  /  'f  n(v(x))dx  +  /  ^  . (Y( v(s) ) )ds,  v  f  w 

n  U  r 

a  continuous  and  convex  function  'f  s  W  ♦  R  with  subgradient 
(u)  *  3*>(u)  «  a0(u)  +  Y  (a_1<Yu)) 

bounded  from  W  to  W*  •  That  is,  F  «  (u)  if  and  only  if  there  exist  f  t  <»0<u)  in 

L2(fl)  and  f  e  a_1(Y(ut)  in  L2(D  for  which 

(7.1)  F(v)  -  /  f  <x)v(x)dx  *  f  t  ( s) v( s) ds ,  v  e  V  , 

a  v  r 

so  the  formal  and  boundary  parts  of  A  are  given,  respectively,  by 

(7.2)  Aq (u)  -  a0(ul>  3a<u)  -  (Yu)  . 

From  Theorem  2  we  obtain  the  existence  of  a  weak  solution  of  the  initial-boundary-value 
problem 


d  2  —1 

yj;  Aq(u)  +  Bq  ( u)  3  fQ  in  L  (0,T;H  ((!)), 

Aou(0)  ,  vQ  , 

S  3a(u)  +  3  (u)  3  g0  in  L2(0,T>h"  ^  (H) , 

[  3Au(0)  3  bQ  . 


This  is  made  precise  in  the  form  (6.5)  and  (6.6)  where  the  operators  are  specified  in 


- - - - - . r-,  . . .  II  . . HI . . . 


Remarks .  By  our  choice  of  V  =  H'(ft),  all  boundary  conditions  in  (7.3)  are  of  variational 
type.  Dirichlet-type  constraints  are  obtained  by  taking  subspaces  of  h'iO). 

We  require  that  fg  and  gQ  be  square-summable  with  values  in  H  '(ft)  and 
H  ^(T),  respectively,  and  we  assume  (6.11)  to  obtain  coercivity  of  B.  The  boundedness 
assumptions  on  (k  =  0,-1)  can  be  relaxed  somewhat  by  using  embedding  theorems,  e.g. , 

of  W  into  LP ( ft ) . 

There  is  no  bound  on  the  degeneracy  permitted  in  the  operator  A ;  we  include  even  the 
(uninteresting)  elliptic  case  A  5  0.  The  case  of  Ag  “  0  leads  to  an  evolution  on  the 
boundary  subject  to  an  elliptic  equation  in  the  interior;  such  problems  arise  from 
diffusion  in  a  medium  bounded  by  material  of  markedly  lower  diffusivity  (25) . 

The  classical  porous-media-equation  and  the  weak  form  of  the  two-phase  Stefan  free¬ 
boundary  problem  are  included  in  (7.3).  In  the  latter,  the  enthalpy  is  given  by 
cigfs)  »  (1  +  cH(s))s  +  LH(s)  where  L  >  0  is  the  latent  heat  of  fusion  and  H(>)  is  the 
Heaviside  function  [14,  16].  Such  problems  arise  in  welding  with  the  nonlinear  term 
3g(u)  representing  a  source  of  heat  due  to  electrical  resistance. 

Note  that  each  solution  of  (5.1)  is  also  a  (spatially  independent)  solution  of  (7.3) 
so  there  is  much  non-uniqueness  in  (7.3). 

(b)  Pseudo-parabolic  Equations.  Here  we  set  V  =  Hj(ft),  so  T  *  {0}  and  all  boundary 
conditions  are  of  Dirichlet  type.  The  operator  A  is  given  as  above  by  (7.2);  the 
operator  8  is  also  given  as  before  but  we  shall  only  assume  (6.7),  not  (6.11).  On  the 
space  V  we  take  the  (equivalent)  scalar  product  and  corresponding  Riesz  map 

n 

Ku(v)  =  f  7  n  u( x)D  v( x) dx,  u,v  t  v  , 

*  **  V  le 

ft  k-1 

n  2  2  *  1  1 

so  we  have  R  =  -i  S  -  7  D.  .  Assume  f  c  I.  (0,T;H  (ft))  and  v„  c  a  (u  ),  u  c  H  (ft) 

n  *-  k  n  oonno 

k«1 

are  given.  Then  from  either  Theorem  l  or  from  Theorem  3  we  obtain  existence  of  a  solution 


of  the  problem 


u  e  H’tO.T/Hpfnn,  u(0 )  -  uQ  , 

V  e  h’(0,T»H  v(O)  -v0  , 

w  <  L2<o,TiH_1<nn  , 

|r  (v(t)  -  4  u(t))  +  w(t)  -  f.(t)  , 
ot  n  o 

v(t)  e  Aq (u( t) ) ,  w(t)  e  Bq (u( t) )  . 

The  operators  A0  and  bq  are  given  by  (7.2)  and  (6.9)  respectively. 

Remarks.  The  partial  differential  equation  in  (7.4)  is  of  the  form  of  a  nonlinear 

3 

parabolic  plus  the  tern  -r—  A  u(x,t).  Such  equations  are  known  to  arise  in  various 

ot  n 

diffusion  problems  and  ere  called  pseudoparabolic  [9,  15,  28].  Similar  problems  with 

variational  boundary  conditions  can  be  considered)  we  obtein  weak  solutions  in  the  form 

• 

(6.4).  However,  since  R  (A.  +  R)  ■  H  (0),  we  cannot  use  Lemma  6.1,  in  general,  to 

g  0 

deduce  (6.6).  This  situation  occurs  even  in  the  linear  case  [26]. 

The  operator  -4  in  (7.4)  can  be  replaced  by  the  Riess  operator  of  any  equivalent 
n 

scalar  product  on  hJ  (!)) .  This  trivial  observation  is  useful  in  introducing  elliptic 
linear  operators  in  its  place. 

we  have  not  made  use  of  the  fact  that  only  one  of  the  operators  A,B  need  be  a 

subgradient.  In  particular,  we  are  free  to  add  to  one  of  A  or  8  any  linear  combination 

of  first  order  derivatives.  (See  Example  (d)  below.) 

Non-uniqueness  of  solutions  of  (7.4)  follows  from  that  of  solutions  of  (5.1). 

In  the  preceding  examples  the  nonlinearity  arises  from  the  local  dependence  on  the 

solution,  e.g.,  from  nonlinear  functions  of  the  values  of  u  or  Vu  at  each  point  of  0. 

we  next  display  examples  of  global  nonlinearity  arising  from  the  "total  energy"  or  the 
"total  flux"  in  the  system.  The  following  preliminary  result  will  be  useful. 

Lemma  7.1.  Let  a(*,*)  and  b(*,*)  be  continuous,  bilinear,  symmetric  and  nonnegative 
real-valued  functions  on  the  Hilbert  space  V.  Then  for  a,B  t  R,  the  function 
l(u)  =  ^max{a(u,u)  ♦  a,b(u,u)  +  6),  u  t  V 
Is  convex,  continuous  and  its  subgradient  is  given  by 

-28- 


r 


(7.4) 


3v>(  u) 


{A(u) } 

t ( XA  +  (l-X)B) (u),  0  <  X  <  1} 
{B( u) } , 


if  a(u,u)  +  a  >  b(u,u)  +  8 

if  a(u,u)  +  a  ■  b(u,u)  +  B 

if  a(u,u)  +  ct  <  b(u,u)  ♦  B 


where  Au(v)  «  a(u,v),  Bu(v)  «  b(u,v),  v  e  V. 

Proof:  We  need  only  to  compute  3tf(u).  For  the  first  and  last  cases  we  compute  the 

Gateaux  derivative  lim  {^<u  +  tv)  -  *{u))/t)  to  obtain  the  desired  results.  Now  assume 
t*0 

a(u,u)  +  a  »  b(u,u)  +  $.  An  easy  computation  gives 

t  '(^(u+tv)  -  <(  u) )  ■  max{a(u,v)  +  a(v,  v)  ,b(u,v)  +  -j  b(v,v)} 
so  we  have  the  equivalence  of  f  e  3^(u), 

f(v)  <  t  ^(^(u+tv)  -  s?(u)>,  v  e  V,  t  >  0  , 


and  of 


f(v)  <  max{a(u,v),b(u,v)},  v  e  V  . 

This  is  equivalent  to  f  =  XAu  +  (1  *  X)Bu  tor  some  X,  0  <  X  <  1. 

(c)  Energy -Dependent  Elliptic-Parabolic  Equation.  We  shall  use  Theorem  2  with  the 

operator  8  given  by  (6.8),  so  we  assume  (6.7)  and  (6.11).  Choose  V“H1(fl)  so  the 

space  of  boundary  values  is  T  »  (T).  Define  on  W  5  L2(£2)  the  function 

*>(u)  “V-jmaxd,  /  |u(x)|2dx},  u  e  w  . 

£2 


The  subgradient  A  *  3v?  is  given  by  Lemma  7.1  and  we  have  A  »=  Ag  ■  Afl,  Tt^( A )  •  L  (0). 
Finally,  let  vQ  e  L2(fl),  fQ  e  L2(0,T:  L2(!?)),  gQ  e  L2(0,T,-h"  ^  O')  >  be  given  and  define 
f(t)(v)  =  /  fQ(x,t)v(x)dx  +  g0(t)(yv),  v  e  V  . 


Then  we  obtain  a  weak  solution  of 


(7.5) 


It  +  B0(u)  *  f 0 

v(x,0 )  ■  vQ(x) 

38(U)  3  g0 


in  L2(0,T»H_1(fl)) 
in  L 2 (0) 

in  L2(0,T;H~  1/2  (T) ) 


where  v  is  determined  by 


if  /  tul  dx  <  1  . 

a 

if  /  |u|2dx  «  1  , 

a 

if  /  |u|2dx  >  1  • 

a 

Thus,  the  type  of  the  equation  is  either  elliptic  (with  parameter  t)  or  parabolic  and 

depends  on  the  total  energy  /  |u|2dx. 

“  1  2 

(d)  A  Flux-Dependent  Equation.  Take  V  =  hJ(Q),  M  ■  L  (S)  and  T  =  {0}.  Let  the  convex 
function  q  and  its  bounded  subgradient  aQ  =  3^  be  given  as  above  in  (a),  and  define 
A  «  a0  in  L  (a);  cf.  (7.2).  Denoting  the  gradient  of  u  by  Vu,  we  define  the 
continuous  convex 

i(i(u)  =  VjmaxtN,  /  |Vu(x)|2dx},  u  e  V  . 

a 

* 

Let  b  f  r"  and  define  8  t  V  *  2V  by 

B(u)  =  b  •  Vu  +  3i|i(u)  . 

Note  that  8  is  maximal  monotone,  bounded  and  coercive.  Let  v0  e  Rg(A)  and 
2  -1 

fQ  e  L  (0,TjH  ( a ) ) .  From  Theorem  2  we  obtain  existence  of  a  solution  of  the  problem 

u  6  L2(0,T;Hg(a)  )  ,  v  e  H1(0,T;H_1(a)),  w  <?  L2(0,T;H~1  (£))  ) 

+  b  •  Vu  -  K( /  |Vu|2dx)i  u  =  fn 
3t  f!  "  0 

v(x,t)  t  aQ(u(x,t),  v(x,0 )  =  Vg(x) 

where  the  maximal  monotone  K  :  R  ♦  R  is  given  by 

f  {0}  ,  S  <  N  , 

I 

KCs)  »  j  (0,1],  S  *  N  , 

i  (1)  ,  s  >  N  . 

Remarks.  In  the  region  where  f  |Vu|2dx  <  N  the  equation  in  (7.6)  is  a  conservation  law 

a 


f 

I 

(7.6)  i 


{0}, 

{Xu  :  0  <  X  <  1} 
{u}  , 


of  the  forn 


(7.7) 


3t  +  b  7g(v)  J  f0 


where  the  maximal  monotone  g  :  R  *  R  is  the  inverse  to  a^.  Thus  (7.6)  suggests  a 
penalty  method  [18]  to  approximate  solutions  of  (7.7).  We  shall  develop  these  observations 
elsewhere. 

In  order  to  consider  (7.6)  in  the  form  (6.1)  it  is  essential  that  8  is  not  required 
to  be  a  subgradient. 

(®)  Elliptic-Parabolic  Systems.  Our  final  example  consists  of  a  pair  of  equations  of  the 
type  given  above  in  (a)  that  are  ( nonlinearly)  coupled.  Pori-  0,1  and  k  «  0,-1,  let 
^  :  8  ♦  R  be  convex  and  continuous  with  subgradient,  a*  =  satisfying 

(7.8)  |w|  <  C< | s |  +  1)  for  w  e  a^(s),  s  e  P  . 

On  the  product  space  W  =  Hr(Sl)  x  Hf((l),  r  <  1,  we  have  the  continuous  trace  operator 

12  12  2  2 
Y([u  ,u  ])  —  [Y(u  )  ,Y(u  )]  which  maps  W  into  L  (T)  x  L  ( I* ) .  Thus  we  define  by 

2  2 

v'(v)  =lf  ( v^  ( x)  )dx  +  l  f  s-‘i.(Yv1(s))ds,  v  -  [v’.v2]  e  w  , 

i=1  S2  i-1  r 


a  continuous  and  convex  function  whose  subgradient  is  given  by 

A(u)  S  3v'  ( u )  «  (a^tu1)  +  Y  (a^  (Ylu1 ) ) ) ,  a2(u2)  +  Y  (a2 1  (Y  (u2 ) ) )  ]  , 

u  -  [u  , u  ]  f  w  . 

«  1  i 

w 

The  operator  A  :  W  -»  2  is  bounded;  its  formal  and  boundary  parts  are  given, 
respectively,  by  (see  (7.2)) 

(7.9)  A0(u)  -  (dgCu1),  a2(u2)],  3^(u>  =  [a^  (Yfu1) )  ,  “^(Yfu2))!  ■ 

Hereafter  we  restrict  Y  to  the  product  space  V  =  H 1  ( VI )  x  H^O).  Assume  we  are  given  a 
set  of  continuous  and  convex  functions  i|j^  s  R  ♦  R  for  i  -  1,2,  -1  <  k  <  n,  whose 
subgradients  3^  =  3'4i*  all  satisfy  the  estimate  (6.7).  For  1  «  1,2  we  define 
^  :  H1 (ft)  ♦  R  as  in  Section  6;  its  subgradient  is  then  given  by  (see  (6.8)) 


6i(ui)  -  SH/Su1) 


=  X  Dksk(Dkui) 

k=0 


+  Y 


3*^ (Yu  ) , 


u  «  H  (0) 
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The  formal  and  boundary  parta  of  8*  are  given  by  (6.9)  and  (6.10)  for  each  of  1  •  1,2. 
Thua  we  have  two  pairs  of  operators  similar  to  the  pair  in  Example  (a).  The  coupling  of 
the  corresponding  equations  will  be  attained  by  a  maximal  monotone  graph  u  :  *  *  2  which 
is  bounded,  i.e.,  (7.8)  holds  for  w  e  U(s).  Then  we  define  a  maximal  monotone  operator 
M  on  R  x  It  by 


Mil.  •  ]>  “  (tw,-wl  :  w  «  y(s1  -  s2)},  ts^l  *  R  *  *  • 

2  2 

This  operator  M  induces  a  corresponding  operator  on  L  (!>}  *  L  (S5>,  hence,  from  V 
into  v*,  which  we  also  denote  by  M.  Finally  we  define 

3([u1,u2J)  »  CiJ 1  ( u^ ),  82(u2>]  ♦  M(u, ,u2>,  lui'u2^  *  v  * 

This  6  is  the  sum  of  maximal  monotone  operators,  each  of  which  is  defined  on  all  of  V, 


.>■ 

jV**-" 


so  6  is  maximal  monotone.  Similarly  B  is  bounded,  and  we  note  that  B  is  coercive  if 
1  2 

both  of  B  and  t>  are  coercive. 

Assume  that  we  are  given  the  following  data: 

fg  t  L2(0,T|H'1(fl)),  L2(0,T»H*1y^  (D),  i-  1,2  , 

[v’.v2]  f  JMA0>,  (v’^vf,)  r  J*giaA>  • 

If  the  functions  {0*  :  -1  <  k  <  n}  satisfy  (6.11)  for  both  i  »  1  and  i  »  2,  then  from 
Theorem  2  it  follows  there  exists  a  weak  solution  of  the  system 

f  ft  *  Bplu’tx.t))  +  u(u\x,t)  -  U2(X,t))  3  fj(x,t)  , 


a2(u2(x,t))  ♦  B2(u2(x,t))  -  u(u1(x,t)  -  U2(x,t)>  3  f 2( x, t)  , 

in  L2(0,T»H_1 (fi))  , 

af1(rui(s,t))  +  3  t(ui(s,t)>  )  gj,  1  ■  ’.2,  in  L2(0,T;H~ 1/j  (T) )  , 
6 


Op(ui(x,0))  3  v*(x)  ,  i  -  1,2,  in  L2(fi)  , 


a^1(Yui(s,0))  »  V^(s)  ,  i  -  1,2,  in  I,2 ( T) 


Remarks.  All  of  the  operator#  In  this  system  are  (possibly)  multi-valued  so  each  of  the 
"equations*  should  be  made  precise  as  was  done  in  our  preceding  examples.  See  (6.9). 
(6.10),  (7.2)  and  (7.9)  for  related  computations. 

The  only  requirement  on  the  is  that  they  be  maximal  monotone  graphs  in  R  which 

satisfy  the  bound  (7.8).  Thus  much  degeneracy  is  possible  in  the  leading  operator 
given  by  (7.9).  Related  stefan-type  free  boundary  problems  can  be  so  considered. 

Interesting  examples  of  the  coupling  term  arise  in  applications  to  diffusion 
problems.  These  include  problems  with  a  semi-permeable  membrane,  u(s)  *  s+  (where  s+ 
denotes  s  if  s  >  0  and  0  otherwise),  or  those  with  a  threshold  phenomena. 

u(s)  -  (s  -  c)+  -  (-s  -  e)+.  The  operator  M  as  given  above  is  a  subgradienti  this  is 
easily  verified  by  showing  it  is  cyclic  monotone  (1].  However  we  may  add  to  M  non- 
symmetric  monotone  terms,  for  example,  [-82.8^1 »  and  thereby  obtain  systems  of  the  form 
(6.1)  in  which  8  is  not  a  subgradient. 

Systems  of  equations  of  pseudoparabolic  type  can  be  resolved  similarly  by  Theorem  1. 
For  example  we  can  choose  V  -  (ft)  x  H^(ft)  with  scalar-product  on  each  factor  as  given 
in  Example  (b)  and  obtain  existence  of  a  solution  of  the  problem 

(a^luVx.t))  -  ^u’lx.t))  ♦  B’lu’tx.t))  +  Mu’lx.t)  -  u2(x,t))  3  fj(x,t), 

(a2(u2(x,t))  -  dnu2(x,t))  +  B2(u2(x,t))  -  u(u1(x,t>  -  u2(x,t))  3  f2(x,t), 

in  L2(0,T!h"1(!1)), 

U-*  f  H '  (  0,T;Hq  (ft) ) ,  u^(x,0)  -  u^tx),  Og(ui(x,0))  3  v^(x),  j  »  1,2, 

in  L 2(ft) 

where  the  data  is  given  as  above  with  v,^  f  AQ(u^)  for  ”  1  • 


j  i 

Vi 
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